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Abstract 

We investigate how the topological charge density in lattice QCD simulations is affected by violations of chiral symmetry 
caused by the fermion action. To this end we compare lattice configurations generated with a number of different actions 
CS| including first configurations generated with exact dynamical overlap quarks. We visualize the topological profiles after 
mild smearing. In the topological charge correlator we measure the size of the positive core, which is known to shrink to 
zero extension in the continuum limit. To leading order we find the core size to scale linearly with the lattice spacing with 
the same coefficient for all actions, even including quenched simulations. In the subleading term the different actions 
vary over a range of about 10%. Our findings suggest that non-chiral lattice actions at current lattice spacings do not 
differ much for observables related to topology, both among themselves and compared to overlap fermions. 
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1. Introduction 

The topological charge density has always been an im- 
portant observable to characterize the vacuum of Quantum 
Chromodynamics (QCD). Instantons are the most natural 
candidates carrying topological charge. Models based on 
them are intuitive, e.g. induce a chiral condensate through 
near zero modes, but also limited, since it became clear 
that the instanton ensemble cannot be dilute in reality. 
Such models are also not complete, as instantons do not 
explain confinement. 

In recent years a more elaborated view has emerged. 
The relevant topological degrees of freedom depend on the 
scale. While instantons are seen in the IR (after strong fil- 
tering), lower dimensional objects are present in the UV, 
for a review see [1 . A laminar structure of them seems 
needed to match the known behavior of the topological 
correlator [2 (see below). This picture is based in partic- 
ular on the analysis of lattice configurations. 

Lattice QCD has seen progress towards smaller lat- 
tice spacings and realistic pion masses. Chiral symmetry, 
however, is realized very differently in these simulations. 
Most of the fermion actions used keep only an approxima- 
tion of it. The reason is that the only known option to 
implement chiral symmetry exactly are Ginsparg- Wilson 
fermions, usually realized as overlap fermions, which are 
numerically very expensive. 



* Corresponding author 
Email address: florian.gruber@physik.uni-regensburg.de 
(Florian Gruber) 



Since chiral symmetry is a key feature of QCD, it is 
not clear a priori, how harmful this approximation is for 
physics results. However, the fundamental aim of lattice 
QCD is to avoid uncontrolled approximations and, there- 
fore, it is necessary to perform benchmark comparisons 
with dynamical overlap fermions to quantify the size of 
systematic errors due to chiral symmetry violation. 

In this Letter we analyze the topological charge in sim- 
ulations obtained with different state-of-the-art lattice fer- 
mion actions. Topology is intimately connected to chiral 
symmetry, both through the U{1) anomaly giving mass to 
the T]' and through the index theorems linking fermionic 
zero modes to topological charge. We, therefore, believe 
that the topological charge density, although it is a purely 
gluonic quantity, is well suited to test the effects due to 
violation of chiral symmetry for different fermion actions. 
The distance at which the topological charge correlator 
changes sign seems to be especially robust against small 
changes in the simulation parameters. 

We use the field-theoretic definition of the lattice topo- 
logical charge and compare topological structures obtained 
after a few sweeps of smearing (alternative filtering meth- 
ods, for instance based on chiral lattice Dirac operators, 
have been shown to largely agree on topological structures 
and to be sensitive of the difference between quenched 
and dynamical configurations, see below). We visualize 
the topological profiles after mild smearing though results 
for different lattice spacing a are difficult to compare. To 
reach more quantitative results, we compare the size of the 
positive core in the topological charge correlator and the 
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amplitude of the contact term. 

For these quantities we compare the results with those 
of configurations generated with overlap fermions. These 
latter simulations are too expensive for comprehensive sim- 
ulations, but they allow for a benchmark test of the de- 
scribed nature on small lattices. 

Our main conclusions are that the differences between 
the actions are not very large, even not for quenched simu- 
lations. The sensitivity to the lattice spacing a far exceeds 
that to the fermion action used. Therefore, observed dif- 
ferences, e.g. between quenched and dynamical fermion 
results should not be over interpreted. Obviously similar 
analyses of other quantities are needed to decide whether 
these conclusions are valid in general. 

2. The topological charge density 

2.1. . . .in the continuum 

The topological charge and its (Euclidean) density are 
defined as 

Q = d^xq{x), q{x) = ^^ e^^p^Tr(^F^^(x)Fp^(x)j . 

(1) 
In this Letter we focus on the two-point function of the 
topological charge density 



Cq{x-y) = {q{x)q{y)), 
which integrates to the topological susceptibility 



Xto 



V 



= j {q{Q)q{x))d' 



X. 



(2) 



(3) 



In the continuum, this correlator is negative for any 
finite distance \x — y\ > 0. This follows from reflexion- 
positivity of the theory and the pseudo-scalar nature of the 
topological charge density [3 . At zero distance, however, 
we obviously have a positive correlator, as Cq{0) = (<7(0)^) 
is nothing but the mean-square topological charge density. 

This property on its own would be strange enough, but 
the topological susceptibility is non-negative and vanishes 
in the chiral limit [3^. Therefore, the positive contact term 
has to (over) compensate the negative contribution of the 
correlator in the spacetime integral in Eq. (Isl). 

The topology of the gauge field nicely reflects itself in 
fermionic observables. Due to the famous index theorems, 
the total topological charge gives the difference of numbers 
of zero modes with left-handed and right-handed chirality, 
Q = ul — ur^ of the chiral (massless) Dirac operator. 



2.2. 



and on the lattice 



The topological charge, as a genuine continuum quan- 
tity, is not defined in a strict sense on a lattice. Of course, 
the expectation is that its properties will be recovered in 
the continuum limit. At finite lattice spacing, however, it 
cannot be defined without ambiguities. 



Throughout this Letter we will use the discretized 
topological charge density of Eq. (IT]) (also called field- 
theoretic charge density) , where F^jj is an improved lattice 
field strength tensor [4 , which combines 1 x 1, 2 x 2 and 
3x3 loops to achieve 0(a^)- improvement at tree-leveQ 
This definition is better behaved in the continuum limit 
than the naive discretization with only 1x1 loops and 
gives topological charges Q = ^^ q{x) closer to integers. 

Calculating the topological charge density on a typi- 
cal Monte Carlo configuration one finds that the signal is 
dominated by short range fluctuations. Thus, an operator 
which couples strongly to these fluctuations is ill-defined. 
Small variations of the gauge links would lead to big vari- 
ations of the resulting topological density. Hence, filtering 
methods are necessary to extract the relevant longer range 
degrees of freedom. 

Many methods have been developed since the advent 
of lattice QCD. Amongst them are cooling, APE smear- 
ing [5^, stout smearing [6 and variants thereof, as well as 
techniques based on the low mode truncation of the lattice 
Laplacian [T. 

A different definition of a lattice topology uses the 
index theorem, which for Ginsparg- Wilson operators is 
valiqj and can be formulated in terms of a local density, 
Qferm{x) = Tr (75 (I)(x, x)/2 — 1)) [H] . This density can be 
filtered naturally by including in a spectral representation 
of D only the lowest IR modes [H |10] (see also [11 for a 
recent filtering method relying on adjoint fermions). 

The structures which emerge after these filtering pro- 
cedures are not unique. They strongly depend on parame- 
ters like numbers of iterations or mode number. However, 
in the weak filtering regime one finds that different meth- 
ods reveal similar local structures if one carefully matches 
their parameters p!2HT5] . 

3. Smearing 

In this work we use the improved stout smearing algo- 
rithm [151 [E] • This is an iterative procedure acting on the 
gauge links. Each step consists of a simultaneous update 
of all links by 



Uf, 



^{x) = ex.p{ipn{x)) • U^{x), 



(4) 



Here p is a parameter which determines the strength of 
smearing and Q{x) is a hermitian traceless matrix built 
from the weighted sum of all 1 x 1 and 2x1 Wilson loops 



■•^For details on the improvement coefficients see '4^ 

^ Still the total topological charge on the lattice is not unique since 

it depends on the solution of the GW relation, e.g. on the kernel of 

the overlap operator. 



fermion (gauge) action 


Nf 


a[fm] 


V[a4] 


m^[MeV] 


Nconf 


Ref. 


twisted mass (Sym) 


2 


0.10 - 0.063 


24^48 - 32^64 


Si 500 


20 


|T71 


twisted mass (Sym) 


2 


0.10-0.051 


20^48 - 32^64 


Si 280 


20 


m 


np imp. clover (plaq.) 


2 


0.11-0.07 


24^48 


Si 500 


20 


m 


np imp. clover (plaq.) 


2 


0.10-0.07 


32^64 


«250 


20 


m 


asqtad staggered (LW) 


2+1 


0.15-0.09 


16^48 - 28^96 


«500 


5 


m 


chirally imp. (LW) 


2 


0.15 


16^32 


Si 500 


20 


m 


top. fixed overlap (Iw) 


2 


0.12 


16^32 


Si 500 


20 


m 


overlap (LW) 


2+1 


0.13-0.12 


12^24 


Si 500 


15 


m 



Table 1: Configurations used in this work. (Abbreviations: Sym = Symanzik , plaq. = plaquette, LW = Liischer-Weisz, Iw = Iwasaki) 



including the old link: 



n{x) = J2 



5-2e 



12 



W^ixi(x) 



W2xl{x) + Wi^2{x) 



hermitean, 



(5) 



where we used |(l^j,(x) - n^{x)) - | Tr{l]],(x) - n^{x)} 
as the hermitian traceless projection. 

While other smearing methods, as well as cooling, de- 
stroy topological structures in the long run, this algorithm 
is designed to preserve instanton-like objects in a wide 
range of their size parameter, if we choose e = —0.25 and 
p = 0.06 [16 . Moreover, 5 steps of this improved stout 
smearing were found to produce a topological charge den- 
sity very similar to the fermionic one, ^ferm [15 , which is 
why we stick to these parameters in this work. 



4. Lattice configurations 

4-1- Overview of lattice configurations used in this work 

We use full dynamical Nf = 2 and Nf = 2 + 1 flavor 
configurations from different fermion formulations. They 
are available through the International Lattice Data Grid 
[23] in a wide range of lattice spacings, lattice volumes and 
pion masses. 

The fermion actions in these simulations may be clas- 
sified with respect to their chiral symmetry. Variants of 
Wilson fermions like twisted mass [17] and nonperturba- 
tive clover [TF break chiral symmetry explicitly (by a lat- 
tice artefact). Staggered fermion actions like asqtad [19 
possess a remnant chiral symmetry. For those actions we 
have enough data to perform an extrapolation a ^ 0. 

We also use chirally improved fermions, which are an 
approximate solution to the Ginsparg- Wilson relation [20 
and as exact solutions we take overlap fermions with an 
extra topology fixing term [21 and full dynamical overlap 
fermions, see next subsection. More details on the con- 
figurations can be found in Tab. [l] The pion masses for 
the configurations used vary between 250 and 500 MeV, 
which does not seem to be a major problem as little mass 
dependence is found. 



4-2. Dynamical overlap configurations 

The overlap Dirac operator [24 at a mass parameter /a 
and for some suitable kernel operator K is 



L)o = (1 + /i) + (1 - /i)75sign(i^). 



(6) 



The principle algorithmic challenge in the generation of 
dynamical overlap gauge field ensembles is the difficulty 
changing topological charge during the Monte Carlo his- 
tory. As in the continuum, to move to different topological 
sectors for a Ginsparg- Wilson lattice gauge theory requires 
a discontinuous change in the action. For an overlap ac- 
tion, where there is a clear and consistent definition of 
the topological index, a topological index change is equiv- 
alent to a small eigenvalue of the kernel operator crossing 
zero. From this viewpoint, there are two ways in which 
topology changes may be suppressed: either the density of 
low kernel eigenvalues is decreased, or there is a potential 
barrier between topological sectors. The Hybrid Monte 
Carlo trajectory can only penetrate this barrier if it has 
sufficient momentum (there is a clear analogy between the 
mechanics of the overlap HMC algorithm and the classi- 
cal mechanics of a particle approaching a finite potential 
barrier: the only difference is that instead of three spatial 
dimensions in classical mechanics, the number of 'spatial' 
dimensions is of order of the lattice volume in the molec- 
ular dynamics). Changing topology becomes harder for 
most lattice discretizations as they approach the contin- 
uum limit. For overlap fermions we have this challenge at 
any lattice spacing. 

JLQCD's topology fixed simulations [21 add extra 
fields to the action which act as a determinant factor 
det{i^^/(i^^ -\- T]^} in the Boltzmann weight. This has 
the effect of suppressing (and removing) the low kernel 
eigenvalues, such that the simulation is carried out in a 
fixed topological sector, which greatly reduces the com- 
puter time required for each trajectory. 

The physical effects of the topology fixing are removed 
in the infinite volume limit [25 , while the suppression 
of low modes is said to mimic the expected behavior in 
the continuum limit. Questions remain, however, about 
the ergodicity of this approach; in particular if the cre- 
ation of pairs of would be zero modes (for example an 
instanton/anti-instanton pair) is driven by mixing between 
low modes of the kernel operator. 



The second possibility is to generate overlap ensembles 
allowing (and even encouraging as much as the underly- 
ing physics will permit) topological charge changes. The 
standard, and only known, method to modify the Hybrid 
Monte Carlo algorithm for a discontinuous action is the 
transmission/reflection algorithm first proposed in [26 and 
refined in several subsequent papers [27H32] , together with 
the eigenvector differentiation algorithm proposed in [33] . 
which is required to accurately differentiate near degener- 
ate pairs of eigenvectors whose eigenvalues lie close to a 
discontinuity. The chief focus of current algorithmic re- 
search has been to improve the rate of topological tunnel- 
ing. It can be shown that if the discontinuity of the action 
at the boundary between topological sectors is A^*, then, 
with an initial Gaussian momentum distribution, the prob- 
ability of tunneling from one topological sector to another 
scales like ~ min(l,e~^'^) for all the possible transmis- 
sion/reflection algorithms. 

To obtain a high rate of topological charge changes and, 
presumably, a small topological autocorrelation time, it is 
crucial to use the form of the action that has the smallest 
A^*. When trying to find it one encounters the problem 
that the pseudo-scalar estimate of the determinant gives 
a very bad estimate of A^ [27 , in particular for small 
fermion mass {AS ~ /i~^). 



y ZAOpseudofermion / 



= /^ti 



1 



^dId,{\>o) dIDo{x<o)^ 



> (A5real> = 2 log (l - 2(1 - m)47^^o), 



Do 



(7) 



where ipo is the kernel eigenvector whose eigenvalue chan- 
ges sign and (j) is the pseudofermion field used to estimate 
the determinant. 

To improve the pseudofermion approach one can intro- 
duce the modified Dirac operator 



L>o = (1 + /i) + (1 - /i)75sign(i^ - Ao), 



(8) 



for a real parameter Aq which remains constant over the 
course of the trajectory. With this operator one can 
rewrite 

det Do = det Do det{Do/Do)- (9) 

The first term det Do is continuous during the topological 
index change and can be estimated using pseudofermions, 
while the second term log det {Do /Do) can be calculated 
without pseudofermions. This factorization leads to a cor- 
rect estimate for the action discontinuity: 

A6' = A log det{Do/Do) = A log det Do = A^'reai- (10) 

This factorization method (together with other, similar, 
approaches) thus gives the optimal approach to tunneling 
through the topological barrier, with the highest transmis- 
sion rate possible while conserving the molecular dynamics 
energy. 



Because deti^o at mass zero satisfies a Ginsparg- 
Wilson relation, and we do not permit the index of det Do 
to change during a trajectory, it is possible to run single 
flavor simulations by using the relation 



det Do ex det (^2 + 2/1^ 



± 



2V2 + 2/i2 



(l±75)sign(i^-Ao)(l±75)), 



(11) 



where the sign is chosen to avoid the zero modes and 
the constant of proportionality is a function of the mass 
and the index of det Do- Because the operator on the right- 
hand side is hermitian and positive definite, it is possible 
to use a standard pseudofermion approach to estimate the 
determinant for a single flavor, without need to employ an 
expensive rational approximation. 

The data used in this study is based on three simula- 
tions on 12^24 lattices using a tadpole improved Liischer- 
Weisz gauge action. Though we can simulate individual 
quarks we choose two degenerate light quark masses and 
the strange quark mass is tuned according to the QCDSF's 
prescription [34 . The pion masses, which are restricted 
by the small lattice volume, were measured to be around 
510 MeV, 560 MeV and 600 MeV, at lattice spacings of 
around 0.12 fm. Over 1000 trajectories were generated for 
the lightest pion mass, and over 600 trajectories were used 
for the heavier pion masses. We used a Wilson kernel with 
one step of improved stout smearing [16 . 

The limiting autocorrelation time for these ensembles 
was, not surprisingly, the topological one. For the lightest 
mass, there were on average 1.05 attempted topological 
index changes per trajectory, of which 37.5% resulted in 
a change of the topological index. This (naively) corre- 
sponds to one topological index change for every 2.8 tra- 
jectories of length 0.5. (This excludes trajectories where 
the topological charge changed twice or where the trajec- 
tory was rejected.) The Monte Carlo acceptance rate was 
around 85%. Values for the other 12^24 runs were similar. 

However, the integrated autocorrelation time for the 
squared topological charge for this lightest mass ensemble 
was around 40 trajectories (of length 0.5), compared to 
around 8 for the plaquette. This problem in global topo- 
logical autocorrelation can be seen in Fig.[l] It can also be 
seen that, although topological index changes occur com- 
monly, there still exist problematic long term correlations 
for the topological index. Obviously this is an issue for 
further algorithmic improvements, which, however, should 
not influence the present analysis significantly. 

5. Results 

5.1. Visualization of the densities 

A direct way of comparing the topological structures is 
to plot three-dimensional slices of the topological charge 
density. It is crucial to notice that a direct comparison of 
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Figure 1: The topological index history for one of the 12^24 dynam- 
ical overlap ensembles. 



the profiles from different actions is only possible if the lat- 
tice spacings are similar because the densities q{x) strongly 
scale with a, namely as a~^. 

In Fig.|2]we show the topological structure of one sam- 
ple configuration for different fermion actions, after 5 steps 
of improved stout smearing. They all have the same lat- 
tice spacing of a = 0.12 fm and the same physical volume 
V = (1.44 fm)^. Below we will discuss how sensitive the 
topological charge density is to a, implying that without 
further analysis one can only compare lattice results for 
different actions if they have the same lattice spacing. For 
comparison we have generated quenched systems with two 
different gauge actions and included their topological pro- 
files in that figure. We have chosen the plaquette gauge ac- 
tion and the Iwasaki gauge action, which are the quenched 
counterparts of the nonperturbative clover action and the 
topology fixed overlap configurations, respectively. While 
this is only a small fraction of the total four-dimensional 
volume, one can see already some of the main properties 
analyzed in detail below. Dynamical lattice simulations 
tend to give larger fluctuations of the topological charge 
density than quenched ones, as already pointed out in, 
e.g., Refs. [I4l|35j. This property is not shared by the dy- 
namical overlap results. Topology fixed overlap fermions 
lie somewhat in between. Below we will quantify these 
observations. 

To stress the importance of the lattice constant a we 
show in Fig. |3] the topological charge density for differ- 
ent lattice spacings but fixed physical volume for twisted 
mass fermions. Two main effects in the continuum limit 
are clearly visible. First, the structure becomes more and 
more fine grained and second the magnitude of the density 
in physical units increases. 

5.2. Topological clusters 

A cluster analysis of the topological charge density is 
a powerful tool to characterize the profile of the structure 
(cf . [121 E] ) • A cluster is defined as a set of neighboring 
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Figure 2: Slices of the topological charge for dynamical overlap 
(ttztt = 600 MeV), asqtad staggered, dynamical overlap with topol- 
ogy fixing term (m-r^ ~ 500 MeV each) and nonperturbative improved 
clover (quite heavy ttZtt ~ IGeV) fermions, all with the same lattice 
spacing a = 0.12 fm. The quenched counterparts of the latter two 
simulation algorithms are also depicted (Iwasaki ^H- top. fixed over- 
lap and plaquette ^H^ clover). The lattice volume is 12^ in all cases. 
This corresponds to a size of approximately (1.5 fm)^ in physical 
units. The color scale is equal in all plots: Blue represents negative 
topological charge and red positive charge. 



lattice points with the same sign of the topological charge 
density. 

To analyze the shape of these clusters, we cut the ab- 
solute value of the topological charge density at a vari- 
able value ^cut and determine the number of clusters 
A^ciuster(<7cut) with \q{x)\ > ^cut as a function of the to- 
tal number of points above ^cut, ^points (<?cut). 

It has been shown that this cluster number obeys a 
power-law and that the exponent 



dlog(7Vcluster(<7cut)) 
dlog(7Vpoints(<7cut)) 



(12) 



of this power-law is highly characteristic for the underly- 
ing topological structure. Pure noise, for example, would 
correspond to ^ = 1 as every point forms its own cluster 
and a very smooth density would have a cluster exponent 
close to zero. 

In Fig. |4] we show this power-law for the actions used 
in the previous section. We find that the power-laws agree 
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Figure 3: Three-dimensional slice of the topological charge density for twisted mass fermions at fixed physical volume V = (1.9 fm)^ after 5 
steps of improved stout smearing. We see an increasing laminar structure for finer lattices. 
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Figure 4: The number of clusters per lattice volume as a function of 
the volume fraction above the cut-off, for the actions used in Fig. |2] 
Fitting a power-law yields an exponent of ~0.85. 



very well and that the exponents are all compatible with 
each other. Therefore, the topological profiles are similar 
and the apparent difference in the three-dimensional vi- 
sualizations in Fig. |2] originates from a difference in the 
absolute values of the densities. 

5.3. Topological charge density correlator 

Another possibility to quantify the topological struc- 
tures is to look at the two-point correlation function of 
the topological charge density. To this end, we compute 
the "all-to-all" correlator of the density (see Eq. ([2|) after 
5 steps of improved stout smearing. At this point we want 
to recall that the resulting density is very similar to the 
fermionic topological charge density [15 and therefore, we 
find the same qualitative behavior for the two-point func- 
tion as in Ref. [36^ . We show an example of the correlator 
for the dynamical overlap ensemble in Fig. [5] 

At small distances the correlator develops a positive 
core of radius Tc, for large distance the correlator is com- 
patible with zero and in between it is slightly negative [36] . 
This behavior is characteristic for all fermion actions and 
quenched ensembles \lb\ . 
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Figure 5: Two-point function of the topological charge density after 5 
improved stout smearing steps for an ensemble of dynamical overlap 
cofigurations. 



5.3.1. Size of positive core 

The radius of the positive core re is given by the (first) 
zero of the correlator. To this end, we interpolate the cor- 
relator between the adjacent data points. Although, there 
are many potential error sources, the position of this zero 
seems to be rather robust and has an excellent signal to 
noise ratio. Most importantly re shows only little depen- 
dence on a when expressed in lattice units, which allows 
to compare configurations generated for different lattice 
spacings. 

We measured this Tc for different lattice spacings, while 
keeping the pion mass approximately constant. Fig. |6] is 
our main result from which we will primarily draw our 
conclusions on the relevance of chiral symmetry for topo- 
logical properties. 

The core size has to vanish in the continuum limit. We 
expect discretization/smearing errors which are indepen- 
dent of the action used and errors related to these actions. 
As all actions are 0(a)-improved the latter should only set 
in at O(a^). Therefore, we fit the data with a function of 
the form 

rc/a = C^B-a^ . (13) 

We find indeed that the coefficients C are for all actions 
compatible with C = 2.15 and that the 5's differ, see 
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Figure 6: Zero of the 2-point function of the topological charge density. The errors - which are partly to small to be visible - result from the 
interpolation between the points of the correlator. We included fits and error bands for twisted mass, staggered and clover action. There are 
three different points from the overlap ensembles at different pion masses: (a) 600 MeV, (b) 560 MeV and (c) 510 MeV. Note that the twisted 
mass results for m-K = 280 MeV and 500 MeV nearly fall on top of one another. 



fermion action 


m^[MeV 


C 


B fm"^ 


twisted mass 


«500 


2.16(2) 


25(3) 


twisted mass 


«280 


2.163(2) 


23.2(8) 


np imp. clover 


«500 


2.148(6) 


3.2(6) 


np imp. clover 


«250 


2.146(-) 


1.6(-) 


asqtad staggered 


«500 


2.14(1) 


12(1) 


quenched Iwasaki 
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2.13(5) 


35(7) 


quenched plaquette 


- 


2.07(4) 


31(4) 



Table 2: Parameters of the single fits. (Only two lattice spacings 
were available for the light nonperturbative clover ensembles and, 
hence, no fit error can be given.) 



Tab. [21 The size of these differences is what we interested 
in. We did not fit higher order terms in the universal part, 
because our data points were not sufficient to do so and 
because we are only interested in the differences between 
the actions. We expected that actions with better chiral 
properties should give results which are markedly closer 
to the dynamical overlap ones than actions with strong 
violation of chiral symmetry. 

Let us further discuss the universality of this scal- 
ing behavior. First of all, if we compare the curves for 
twisted mass and clover fermions at m^ ^ 500 MeV and 
m^ ^ 250 MeV, we find a relatively weak dependence on 
the mass. 

Secondly, we have found the constant C to depend on 
smearing (not shown). This is not surprising as the posi- 
tive part of the topological charge correlation gets smeared 
out over more lattice spacings. 

Fig. |6] also includes the results for dynamical overlap 
fermions with and without topology fixing term at similar 



lattice spacings. The data points for dynamical overlap 
configurations are not on one line, because they belong to 
different pion masses and, hence, different fit curves. 

We have also included two quenched results for com- 
parison. The extrapolations yield a consistent value for the 
constant C with a slightly larger error than the dynami- 
cal configurations. The quenched results are in accordance 
with the results of Horvath et al. [36 . They used the full 
fermionic definition of the topological charge density (not 
truncated in Dirac modes) and find a core size Vc ~ 2a for 
the Iwasaki gauge action. 

The non-universal term Bo? measures how Vc/cl from 
different actions converges towards C. The spread at cur- 
rent lattice spacings a < 0.15 fm is on the order of 10%. 
The overlap data lie in the middle of that range and even 
the quenched results are not far off. Our conclusion is that 
the different actions with their different treatment of chiral 
symmetry do not differ much with respect to this impor- 
tant topological observable. Consequently, all of them can 
be used to study topological quantities. To the extent that 
topological properties are relevant for hadron properties 
this observation also ties in with the fact that quite often 
the differences between quenched and dynamical results 
for these are not huge (typically of the order 10-30%). 

5.3.2. Contact term 

As explained in the introduction, the two-point func- 
tion has to develop a positive contact term in the contin- 
uum limit. Therefore, the mean-square value of the topo- 
logical charge density (^(0)^) has to be divergent in the 
continuum limit. Fig. [7| shows a double logarithmic plot 
of (<7(0)^) in physical units versus the lattice spacing. The 
linear behavior indicates a power-like divergence for a ^ 
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Figure 7: Maxima of the two point function of the topological 
charge density. Dynamical overlap results are labeled by ttZtt as: (a) 
600 MeV, (b) 560 MeV and (c) 510 MeV. Note that the twisted mass 
results for ttZtt = 280 MeV and 500 MeV fall on top of one another. 



for all actions: 



{q{Of) (X a" 



(14) 



for some positive number c. 

This exponent is similar for the different actions and 
its value is around —6. Only the nonperturbative clover 
action deviates from this value with an exponent around 
—7. As the contact term is highly divergent in the con- 
tinuum limit, we do not expect the fitted values to agree 
and, therefore, we do not present them. 

The contact term of the dynamical ensembles is big- 
ger than for their quenched counterparts, but it also im- 
portant to notice, that the contact term for dynamical 
overlap configurations at a = 0.12 fm lies between the two 
quenched simulations (plaquette and Iwasaki). Thus we 
can conclude, that dynamical fermions generate indeed 
larger fluctuations in the vacuum at finite lattice spacing, 
as argued in [S^ , but also that the differences between dif- 
ferent actions are as large as those between quenched and 
dynamical simulations. 

6. Summary 

We have investigated the topological charge density for 
state-of-the-art lattice actions with dynamical fermions, 
including new dynamical overlap simulations. This quan- 
tity was chosen because of the intimate connection between 
topology and chiral symmetry. Different fermion actions 
do generate different topological landscapes, as visualized 
in Fig.[2]and quantified through the topological charge cor- 
relator. The change in the topological observables is not 
very large. The radius of the positive core of the topolog- 
ical correlator Vc approaches zero with the same slope C 
for all actions. In the next-to-leading order re differs, but 
the spread is below 10% and even quenched simulations 



do not produce markedly different results. In particular, 
simulations with exact overlap fermions give results which 
are quite similar to those obtained with topology fixed 
overlap fermions. The differences between quenched and 
dynamical simulations are not larger than those between 
different dynamical fermion actions. Also, the topologi- 
cal charge density seems to be little affected by changes 
in pion mass. In contrast, the effects for (<7(0)^) are large 
but unsystematic. These results are very sensitive to the 
lattice spacing a, implying that one should be very careful 
not to jump to conclusions when comparing topological 
properties of different configurations. If we use our dy- 
namical overlap results as benchmark for the quality of 
the other actions with respect to chirality we have to con- 
clude that all of them are reasonable successful and none 
of them seems to be clearly superior. The differences be- 
tween results for dynamical overlap fermions and topology 
fixed overlap fermions are especially small, as one might 
have expected. 
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